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, The radiative corrections to elastic electron-proton scattering are ana- 

lyzed in a hadronic model including the finite size of the nucleon. For initial 
,-0 ■ electron energies above 8 GeV and large scattering angles, the proton vertex 

ry ^ correction in this model increases by at least two percent the overall factor by 

which the one-photon exchange (Rosenbluth) cross section must be multiplied. 
CSJ ! The contribution of soft photon emission is calculated exactly. Comparison 

is made with the generally used expressions previously obtained by Mo and 
J> ■ Tsai. Results are presented for some kinematics at high momentum transfer. 
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§ ! I. INTRODUCTION 



Electron scattering at intermediate and high energies has been one of the most useful 
means of investigating nuclear structure for over forty years. With the advent of CW accel- 
erators and high resolution detectors such as MAMI and TJNAF it has become clear that 
one must have an accurate estimate of the radiative corrections if meaningful cross sections 
are to be obtained from the experimental measurements. Depending on the experimental 
^ ■ conditions - initial beam energy, momentum transfer, and detector resolution or missing 
mass for the observed particles - the radiative corrections can be as large as 30% of the 
uncorrected cross section. To obtain cross sections which are accurate to 1%, one must then 
know the radiative correction to 3%. 

The theoretical expression for the radiative correction which has been used in the analy- 
sis of almost all single arm elastic electron scattering experiments with beam energies below 
approximately 25 GeV (for which W and Z exchange are in general not significant) is that 
given originally by Tsai [H, in connection with experiments at Stanford, SLAC and 
CEA. That work involved approximations that were both purely mathematical (made in 
performing the integrations needed to evaluate the inelastic cross section) and approxima- 
tions denoted here as "soft-photon approximations" that are directly related to the physics 
in that the effect of proton structure was neglected; in considering the proton legs, only the 
soft virtual (infrared) photon contribution is calculated exactly - approximations are made 
in the hard virtual photon (non-infrared) contribution. In particular, the proton structure is 
neglected by setting the photon momentum square k 2 = in the proton form factor, F(k 2 ), 
thus simplifying the calculation considerably. 
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The purpose of the present paper is to study the radiative correction to elastic electron- 
proton scattering including the internal structure of the nucleon. For this we have considered 
a simple model in which the proton current is taken to have the usual on-shell form. The 
model dependence of the radiative correction is clearly an important question for the analysis 
of electron scattering experiments at the 1% level. This work is an initial study to examine 
the size of internal structure effects. 

The present calculation differs from that of Tsai ffl, in three substantive aspects. 
First, we evaluate the inelastic cross section (emission of soft real photons) without any 
approximation; the relevant integrals have been given in closed form by t'Hooft and Veltman 
H]. In fact, these exact expressions are simpler in form than the approximate ones given in |I[ 
and f2|. We note in particular that in the limit of the target mass M — > oo, corresponding 
to a static Coulomb potential, we obtain exactly the result first given by Schwinger ||. 
Second, in the evaluation of the contribution of the box and crossed box diagrams to the 
elastic cross section we make a less drastic approximation than that made in |T[|. Specifically, 
in the integrands corresponding to the relevant matrix elements, M 2 and M 3 (Eqs. ( |A5|) and 
( |A6|) ), we make a soft photon approximation (setting k = or k = q) in the numerator (as in 
U), but not in the denominators. Again, the required integrals (scalar four-point functions) 
have been given in [EJ; the resulting expressions are again considerably simpler than those 
obtained in ]I] , where the soft-photon approximation is also made in the denominators of M 2 
and M 3 . Finally, in evaluating the proton vertex correction, we have made no soft photon 
approximation for the virtual photon (as was done in [|lj) and have included form factors 
for the proton, taking the proton current to be that indicated below in (|2.4| ). 

The organization of the paper is as follows: In Sec. II we discuss questions concerning 
the electromagnetic nuclear current operator used in this analysis. In Sec. Ill we give details 
of the calculation of the matrix elements and cross section for elastic scattering, retaining 
terms of order a relative to the Rosenbluth (one photon exchange) cross section for elastic 
scattering. Integrals needed for the evaluation of the various matrix elements are written 
explicitly and expressed in closed form in terms of Spence functions (dilogarithms). Details 
are given in the Appendices. 

In Sec. IV we consider the inelastic cross section in detail; as with the elastic cross 
section given in Sec. Ill, the result is expressed in closed form in terms of Spence functions. 
In Sec. V we add the elastic and inelastic cross sections, giving both an analytic expression 
and a numerical evaluation of the radiative correction for various values of the pertinent 
parameters, (initial beam energy, final electron detector resolution, and target nucleus). We 
compare the values of the radiative correction calculated here with those given in M and 



II. ELECTROMAGNETIC NUCLEON CURRENT OPERATOR 

We essentially follow in this paper the convention of Bjorken and Drell ||. The metric 
used is defined by 

Pi ■ Pj = EiEj - Pi ■ Pj (2.1) 

Further, a = e 2 /47r = 1/137.036; m is the electron rest mass; M is the target nucleus 
rest mass; Z the charge of the target nucleus; k the anomalous magnetic moment of the 
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proton; p\ and p^ the initial and final electron four-momenta, respectively; P2 and P4 the 
initial and final target nucleus four-momenta, respectively; q = p\ — p^ = p^ — P2 is the 
four-momentum transfer to the target nucleus for elastic scattering. It will prove useful to 
define, in addition, 

P=P4+P2, Pm=Pl+P3, (2.2) 

from which p 2 = —q 2 + 4M 2 and p 2 m = —q 2 + 4m 2 . Further, we define 
x = (p + Pl ) / (p - Pl ) = (p + Pl ) 2 / 4M 2 

X m = {Pm + Pl) / (Pm ~ Pl) = (Pm + Pl) 2 / 4m 2 (2.3) 

with p\ = —q 2 . In the lab system we have: p\ = (ei,pi); p% = (63, P3); P2 = (M, 0); p^ = 
(M + w,q); u = -q 2 /2M. 

With the aim of presenting expressions which correspond to the experimental conditions 
of high energy electron scattering, we neglect, in the final expressions given in this paper, 
terms of relative orders m 2 /e 2 , m 2 /(— q 2 ), and m 2 /M 2 . Neglect of these terms defines our 
high energy approximation. No assumption is made, however, with regard to the magnitudes 
of M/e x , M/e 3 , or M 2 /(-q 2 ). 

At low momentum transfer the internal structure of the nucleon can safely be neglected in 
the determination of the radiative corrections in electron-nucleus scattering. However, with 
increasing energies and momenta this is in general no longer the case. One of the objectives 
of this paper is to investigate this in a model for the e.m. interaction of a non-pointlike 
nucleon. The most general e.m. off-shell nucleon vertex can be characterized by 6 invariant 
functions As the most simple model we may consider a vector dominance-like model 

for the nucleon current, characterized by only two form factors which depend only on the 
four-momentum square of the photon. It is given by 

T^Ftfhr + KFitf)^, (2.4) 
where the form factors F\(q 2 ) and F 2 {q 2 ) are taken to have a monopole or dipole form: 

F^q 2 ) = F 2 (q 2 ) = C^^j , n=lor2 (2.5) 

with A being a constant of the order of 1 GeV/c. Furthermore, q = p' — p, p and p' being 
the momentum of the initial and final nucleon. Although the quantitative predictions of 
the radiative corrections are expected in general to be dependent on the details of the 
nucleon model assumed, one should already be able to see most of the salient features in 
the present model study. In particular, identifying regions in phase space where the finite 
size of the nucleon may play an important role in the size of radiative corrections can be 
important. In this way one may hope to get some feeling on the reliability of neglecting 
the internal structure of the nucleon as is usually done. The present study is intended as a 
first exploration of the sensitivity on the non-pointlike nature of the e.m. hadronic current. 
As in [0J, although we are primarily interested in electron-proton scattering, the radiative 
corrections studied here can also be applied to electron-nucleus scattering, with appropriate 
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changes in Fi, F2, k, and M. However, even in the case of electron-proton scattering, the 
factor Z is convenient for identifying the contributions from the various diagrams. 

It should be noted that the dressed vertex function, A^, with Eq. (|2.4|) as e.m. current 
operator containing the form factors F n , satisfies a Ward-Takahashi identity 

q^ = F 1 (q 2 )[S- l (p')-S- 1 (p)] (2.6) 



where S is the dressed nucleon propagator. As a direct consequence of (|2.6|), one gets for 
on-mass-shell nucleons, the current conservation 

< p'|A> >= 0. (2.7) 

Obviously, the radiative corrections will in general be sensitive to the choice of the e.m. 
current. Although interesting in its own right we will not address in this paper the issue of 
the dependence of the predictions on this ambiguity. 

In the study of radiative corrections we may distinguish between the elastic and inelastic 
contributions, the latter being the real soft photon emission processes from both the electron 
and hadron. The elastic electron cross section can be determined immediately from the total 
scattering amplitude A4 through the well-known expression 

da= 7( ™f =E / l-M| 2 (2vr)V(p 4 + P3 - P2 -Pi) 



spins 



md 3 p 3 Md?p A 
X (2^ (2vr)3 e4 ' V*> 

For single-arm experiments with unpolarized electrons in which the final proton is not ob- 
served, da must be averaged over initial spins, summed over final spins, and integrated over 
the final proton four-momentum. Up to order a 2 we have for the total scattering amplitude 

6 

M = J2 M m (2-9) 

71=1 

where the various terms correspond to the Feynman graph contributions shown in Fig. [I]. 
Mi is the matrix element for the one-photon exchange diagram 

:«) 

q l + ie 

Its square gives, for high energy electrons, the Rosenbluth cross section 



Mj = Ze 2 u(p 3 ) lfM u( Pl )^-u(p 4 )T^q 2 )u(p 2 ). (2.10) 



da a2 cos2 f 



(*? -|S^)-it(^ + ^ 2 ) 2 tan 2 | 



dtt ' 4e^sin 4 (^/2) ' 

where r] is the lab system recoil factor: For ei >> m, 63 >> m, i] = ei/e% = l+(ei/M)(l — 
cos#) with 9 being the electron scattering angle. We note, in particular, that 1 < rj < x. 
Furthermore, M2 and M3 are the matrix elements for the box and crossed box (two-photon 
exchange) diagrams. M4 is the vacuum polarization diagram (only an electron-positron 
loop is indicated in the figure, but the contribution from higher mass lepton loops can be 
included without difficulty- see ( |A8| ) -( |A12| )). M 5 is the electron vertex correction, and 
M 6 is the proton vertex correction. For completeness, we list in Appendix A the explicit 
expressions for the various Feynman diagrams shown in Figs. and |2|. 
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III. ELASTIC CROSS SECTION 



To evaluate the various one-loop corrections to Eq. ( |2.9|) some tedious algebra has to be 
carried out. We outline the procedure used to evaluate the matrix elements needed for the 
radiative correction to the elastic cross section, M 2 through M 6 , (Eqs. (|Al| ) - (|A10|) ). 



A. Proton vertex correction 



We begin with the matrix element for the proton vertex correction, M 6 , given by Eqs. 
and ( |A4|) ; the much simpler matrix element for the electron vertex correction, M 5 , 
(Eq. (|A1|)), can be deduced quite easily from that. In (|A"4T), each of the three I~"s, given 
by Eq. ( |2.4| ), contains a term with 7^ (which we denote by g) and a term with a^ v (which 
we denote by s). The proton vertex correction A^(p4,p 2 ) then consists of eight terms, 
which we represent symbolically by ggg, gsg, gss, etc. As may be seen after rationalizing 
the propagators, the k dependence of the numerators for ggg, gsg, ... is such that there are 
at most four factors of the form Moreover, the terms with three or four factors ^ may, 
with only a minimum of algebra, be written so that two of these factors are adjacent, giving 
1/tJ/t = k 2 . Although the calculation can equally well be carried out with F\ and F 2 distinct 
functions, we assume F\ — F 2 — F, which simplifies the algebra. The terms ggg, gsg, ... can 
then be expressed in terms of the integrals 



where 



d A k 



- F\k 2 ){l; k^, k,k v ; k 2 ; k,k 2 ; k^k u k 2 ; (k 2 ) 2 }/D(X 2 ) 



(3.1) 



D(\ 2 ) = (k 2 - A 2 + ie)(k 2 -2k-p 2 + ie){k 2 -2k- Pi + ie) 



(3.2) 



For form factors having the form given in (2.5), the integrals in ( p.l|) could all be evaluated 
as indicated for three-point functions in 0, Sec. 5, and ||, Appendix E. However, in the 
interest of obtaining a relatively compact analytic expression in closed form, we have used 
an alternative procedure. As given here in Appendix B, the integrals may be expressed 
in terms of their moments, defined by Eqs. ( |B4| )-([B"6D and ( [B10| ). After straightforward 
though somewhat tedious algebra, the terms ggg, gsg, ... are then expressed in terms of these 
moments. Next, for form factors of the form given in fl2.5| ), we show that all of the moments 
may be expressed in terms of three functions, <pk, which obey a three-term inhomogenous 
recursion, and this is used for their evaluation. Finally, we note from ( |B51| ) - ( |B57| ) that the 
terms ggg, gsg, ... may be usefully grouped by writing them in the form 



(9 + s)g(g + s) = F{q 2 



Gx(q 2 h, + G 2 (q 2 ) 



2M 



(3.3) 



and 



(g + s)s(g + s) = nF(q 2 



A 1 (g 2 ) 7M + A 2 (g 



2n «V9 



2M 



(3.4) 
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We note in the expressions for ggg, gsg } ... in Appendix B that the infrared divergent 
terms are all contained solely within ggg and gsg. These are the terms with a factor 0i(A 2 ) 
in ( |B51| ) and ( |B52| ). Since these are precisely the terms which are retained in the proton 
vertex correction in [p]], we separate them for the purpose of comparison with that work, 
writing Mq in the form 

M 6 = M 6 (0) + M 6 (1) (3.5) 

where 

(0) = _^! (2M 2 _ g 2 )0i(A 2 )Mi (3 g) 

The function 0i(A 2 ), defined by ( |BT7|) , (|E25| ) and (|B26|) and evaluated in ( |B33| ) and (|E3T| ), 
is simply related to the function K(p 2 ,p4,) defined in |J by 

K(v p ) = t — 

W*>P3) _ m 2 J (k?-\* + ie)(k 2 -2k-p i + ie)(k 2 -2k-p j + ie) 

viz., 

K(P2,P4) = 2P2 -P40i(A 2 ) (3.7) 



B. Proton self energy correction 



We next consider the contribution of the proton self energy diagrams. It is given by E 
where 



-Tr 



dp 1 



(3.8) 



in which S is the lowest order self-energy contribution 

d 4 k 1 1 



(3.9) 



(2tt) 4 k 2 - A 2 + ie v ' - # - M + ie 
Using the Ward-Takahashi identity Eq. ( ^.6[ ) we find that 

£' = G l (0) 

where G\(q 2 ) is the coefficient of 7 M in ggg ( |B51|) . Explicitly, Gi(0) is given by Eqs. (|B29|) 
- ( |B32|) and (|B38| ). The addition of this contribution to the lowest order vertex correction 
modifies the expressions for (g + s)g(g + s) and (g + s)s(g + s) given above in ( |3.3| ) and 
so that we now have 



(9 + s)g(g + s) = F(q 2 



(G 1 (g 2 )-G 1 (0)) 7 , + G 2 (g 2 )^ 



(3.10) 



and 
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(g + s)s{g + s) = kF(q 2 



X 1 (q 2 )^+(X 2 (q 2 )-G 1 (0)) 



2M 

We then have, for the matrix element including self energy diagrams, 



(3.11) 



M 6 = Nlf + Mf (3.12) 



where 



M ( 0) = o& [ _ K ^ p j +K ^ p j ]Mi (3 13) 

ZlT 

which is the expression given in U, eq. (11.12). The infrared divergent part of these terms 
is cancelled exactly by the infrared divergent terms in the inelastic cross section. The 
contribution of the matrix element which depends on the proton form factor, will be 

considered after we write the electron vertex and box diagram corrections. 



C. Electron vertex correction 

The electron vertex correction, M 5 , may be obtained directly from the proton vertex 
correction, M 6 . The expression A M (p 3 ,pi), Eq. (|A3|), follows from A M (p 4 ,p 2 ) if we retain 
only the term ggg, set F = 1, replace p 2 , Pa and M by pi, p 3 and m, and, after performing 
the integrations, take the limit A — > oo (note Eq. fl2.5j )). Note that p = p± + p 2 is then 
replaced by p m = p 3 +p±, and x, defined in section II, is replaced by x m . Details are given 
in Appendix C. We find 

999 = Gt\q 2 h, + G 2 e \q 2 ) l -^f (3.14) 



where 



2tt [ V 2p m pi J 4 2 \m 2 



G (e) {q 2 )= *i^ lnx \ (316) 

2vr [p m pi J 
Adding the contribution of the electron self energy diagrams gives 

999 = (cfV) - Gi e) (0)) 7, + G 2 e \q 2 )^f (3.17) 

from which 

a f/ r . . T ., . 3p? + 8m 2 , \ / 2m 2 , \ i^g* 
^= — \[ -K(p l ,p 3 ) + K(p 1 ,p 1 ) + -^ lnx m -2 7/ ,+ \nx) 

(3.18) 

For large momentum transfers, — q 2 ^> m 2 , this reduces to 
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999=^ (Pi,p 3 ) + K(p uPl ) + ^ In (-jLj - 2| 7m (3.19) 

Comparing (|2.10|) and ( |AT|) , we then have, for — q 2 ^> m 2 , 

M 5 = ^ |-^bi,p 3 ) + K(p 1 ,p 1 ) + -In f - 2} Mi (3.20) 

which is the expression given in |l[], eq.(II.5). We note that the infrared divergence is 
contained entirely within the terms —K(pi,pz) + K(jpi,pi). The infrared divergent part of 
these terms is cancelled exactly by the infrared divergent terms in the inelastic cross section. 



D. Box and crossed-box diagrams 



The matrix elements for the box and crossed-box diagrams, M 2 and M 3 , are given in 



and (|A6|) . After rationalizing the propagators, the required integrals can, for form 
factors of the form ( |2.5| ), all be written in terms of four-point functions; in principle they 
can be evaluated using ||, Sec. 6, and ||J], Appendix E. For the present work, however, we 
have chosen to evaluate these matrix elements in an approximate manner, but one which is 
less drastic than that employed in We note first in M 2 and M3 that the integrands in 
M 2 and M3 have two infrared divergent factors, [(k 2 — A 2 + ie)((k — q) 2 — A 2 + ie)}" 1 . The 
integrands are thus peaked when either of the two exchanged photons is soft, and become 
divergent when k — > or when k — > q. We therefore evaluate the numerators in M 2 and M 3 
at these two points but make no changes to the denominators. A simple calculation shows 
that in fact that each of the numerators has the same value for k = as for k = q, viz., 
Aipi ■ P2 q 2 M\ in the case of M 2 and 4zp3 ■ p 2 q 2 M\ in the case of M3. We then take this 
factor outside of the integral and are left with a scalar four-point function to evaluate. The 
result has been given in ||, Sec. 6 and Appendix E (b) and is expressed simply in terms of 
logarithms: 

M 2 = -^-Sl >n ( In (4-) Mi (3 21) 

and 



71 IpjJ \ m J \ A 



M 3 = 2? * In ( e -l±M) h ( Z|!) Ml (3.22) 
n |Psl V m J \ A 2 J 

By contrast, in |IJ, in addition to the approximation just described, a soft-photon ap- 
proximation is made in the infrared denominators: Specifically, when k = the factor 
{k — q) 2 — A 2 is set equal to q 2 — A 2 and when k = q the factor k 2 — A 2 is set equal to q 2 — A 2 , 
thus giving two terms and reducing the four-point function to three-point functions: 

* = -|Vfe, (3-23) 

and 
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M 3 = — [K(p 2 ,p 3 ) + K(p 4 , Pl )} M 1 
In 



(3.24) 



(Note eqs.(II.9) and (11.11)). The infrared divergent terms (those with a factor In A 2 ) are, 
for M 2 , the same in ( ggjp and (ggg) , and, for M 3 , the same in (ggp and ( ggg ) . However, 
( p.23| ) and ( |3.24| ) differ significantly from ( |3.21| ) and ( |3.22| ). These latter expressions are 
functions of the momentum transfer, q 2 . The integrals K(pi,pj), on the other hand, are 
functions only of the scalar invariants p 2 , p 2 and Pi ■ pj. In (3.21) and (3.22), M 2 and M 3 
therefore depend only on the initial and final electron energies, and not on the momentum 
transfer ( p 2 ■ Pl = p 4 ■ p 3 = e x M; p 2 ■ p 3 = Pi • Pi = e 3 M). 

From (pD-flMOl), ( gTg ), ( gTg) , (|3^0|) , ( gg| ) and ffT23 ), we can now write the square 
of the matrix element for elastic scattering, including the radiative correction to order a. 
Assuming — q 2 ^> m 2 , and including only electron-positron pairs in the vacuum polarization 
matrix element, we have 



\M\ 2 = | Mi | 



1 + a [f In (^) - f - Kfaps) + K(p 1 , Pl 



_2^ lnr?ln( _|: 



+ ^f[-K(p 2 ,p,)+K(p 2 ,p 2 )} 



A 2 



(3.25) 



+2Re {MjM^ 



E. Contribution of proton form factor 



Finally, we consider the contribution of the term 2Re jM/Mg^j, coming from the inclu- 
sion of form factors for the proton and integration over the entire range of four-momenta 
of the virtual photon in the proton vertex correction. Equations (|3.3p and (|3.4j) define the 



functions Gi(g 2 ), G 2 (q 2 ), X^q 2 ), and X 2 (q 2 ). From ([pp we may write M 6 (1) = M 6 - M ( 



r(0) 



i.e., the term Mq is obtained from the full proton vertex correction by subtracting the 
infrared divergent matrix element which is independent of the proton form factor. We 
therefore define G 1 (g 2 ) and X 2 (q 2 ) to be the expressions G\{q 2 ) and X 2 (q 2 ) from which we 
have omitted the terms with factor 0i(A 2 ) We then write, from (|3.3|) and 



(g + s)g(g + s)' = F(q 1 



(G' l {q 2 )-G , M) 1 , + G 2 (q 



2 ^OjmA£_ 

2M 



(3.26) 



and 



(g + s)s(g + s)' = nF(q 2 



X 1 (q 2 ) 1 ,+ (x' 2 (q 2 )-G' 1 (0) 



2M 



(3.27) 



We now define the functions Fi g (q 2 ) , F 2g (q 2 ) , Fi s (q 2 ) , and F 2s (q 2 ) by 



F(q 2 



G , 1 {q 2 )-G' 1 {Q)) 1 » + G 2 {q 2 ) 



2M 



Fi g (q 2 )^ + KF 2g (q 



2M 



(3.28) 
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*F(g 2 ) [X^h, + {x' 2 {q 2 } - G[(0)) ^ 

Then with the further definitions 

Fi(q 2 ) = F lg {q 2 ) + F ls (q 2 ) 

F 2 (q 2 ) = F 2g (q 2 ) + F 2s (q 2 ) 



F ls (q 2 h, + ^F 2s (q 2 ) 1 -^ 



T^ = F 1 (q 2 ) llM + KF 2 (q 2 



2M 



we have (apart from factors) 



= ^ <p 3 |7>i> <P4|I> 2 



and 



2Re <{ M/iVf^ 



aZ 2 



7T 



2tt 



((psIVIpi) (P4|r,|p 2 )) t ( P3 |7 M bi) (P4|r>: 



This has the same form as 

m\m x = (( P3 \y\pi) (P4|r,b 2 )) f ((paiybi) (P4|r M |P2» 



(3.29) 

(3.30) 
(3.31) 

(3.32) 

(3.33) 

(3.34) 
(3.35) 



with the exception of the replacement T n — > in the right-hand term. Thus, in place of 
the Rosenbluth cross section, (|2.11| ), obtained from J2 spins MlM 1 , we have 



2Re [mJ 



a 2 cos 2 | 



/ spins 

'aZ 2 



spins 



4e 2 7/sin 4 (#/2) \ tt 



where 



FF- 



K 2 g 2 



F 2 F 2 I — 



Q 



(F 1 + kF 2 ) (Fi + KFajtan 







(3.36) 



(3.37) 



4M 2 7 2M 2 

The purely elastic cross section, including radiative corrections to order a, can thus be 
written as 



(3.38) 



where 

Jf(0) 



+ 



and 



a 

7T 

2aZ 

7T 

aZ 2 

7T 



?(!) 



hi 

In 77 In 



-3' 



/ 2 



13 



£4 
ip 4 i 



In x — 1 In 



V A 2 



+ 



9 2 



7T 



(3.39) 



M 2 



7T 



lnxln + -ln'x + 2L( ) + — 



x 



aZ 2 I 



^tFoF> 



£1 

2.W 



(Fi + kF 2 ) ( + kF 2 ) tan 2 f 



7T 



F? 



1 771: 

4Af 2 " r 2 



g 2 
2M 2 



-(F + /tF 2 ) 2 tan 2 | 



6 



(3.40) 
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IV. INELASTIC CROSS SECTION 



In this section we calculate the inelastic cross section, i.e., the contribution of soft photon 
emission by the initial and final electron and proton to the radiative correction. The relevant 
diagrams, with corresponding matrix elements M&i and M b2 , are shown in Fig. |2]. These 
matrix elements are given by 



M, 



61 



iZe (2n) S (p 3 + Pa + k - p\- p 2 ) 



mM 



xu(p 3 ) 



$3 + ft — m + ie 
xu(p 4 )T^u(p 2 ) ' 



^J2uje l e 3 e 2 e 4 
1 

fx — ft — m + ie 



u(pi 



(pi -P3 - k) 2 + ie 



(4.1) 



M b2 = iZ 2 e 3 (27r) 4 8 4 (p 3 + P 4 + k- Pl -p 2 ) 

1 



mM 



xu(p 4 



$4 + ft — m + ie 



v /2uje 1 e 3 e 2 e 4 
1 



-J 



u(p 2 )- 



j) 2 - ft - m + ie J (pi — p 3 ) 2 + ie 



(4.2) 



Making the soft photon approximation, we rationalize the denominators and drop terms 
of relative order k in the numerator and denominator (but not in the delta function), giving 



M M + M b2 = -iZe 6 (2Tty5\p 3 + p 4 + k- Pl -p 2 



mM 



^J2uje 1 e 3 e 2 e 4 q 2 



s / \ — / \p(i / \ / P3 ' ^ Pi-e p 4 -e p 2 -e 

xu(P3n»u(Pi) u(p 4 )T^u(p 2 ) - Z + Z - 

p 3 ■ k pi ■ k Pa - k p 2 ■ k 



(4.3) 



The cross section for soft bremsstrahlung then follows by squaring the matrix element 
Mf,i + Mf, 2 , dividing by the incident flux and the transition rate and multiplying by the 
number of final states. Summing over photon polarizations, we then have 

Z 2 e 6 m 2 M 2 ^ f d 3 p 3 d 3 p 4 d 3 k 4 4 
d(J b = ~ / 9 v 9 n— , 2 f=f 2^ / ~ ~^7T^ 6 (P3 + Pi + k - Pi - p 2 ) 

x \u(P3hAPiMp^u(p 2 )\ 2 j\- 2 (-**- ~~~~T + Z ^h) 2 I 4 - 4 ) 

{Q ) \P3 ■ k pi-k p 4 -k p 2 -k J 

The range of integration in the above expression is determined by the experimental 
conditions. We assume, as in |J, that the final proton and emitted photon are undetected; 
the range of integration in energy and angle of the final electron is determined by the 
entrance slit and spectrometer. We integrate first over d 3 p 4 and are then left with a single 
delta function relating the variables of k and P3: Writing 



d 3 p 4 
~2e~7 



J de 4 5(pl - M 2 ) d 3 p 4 = J d 4 p 4 5(pl - M 2 ) 0(p° 4 ) 



(4.5) 
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with 



t = p 1 +p 2 -p 3 =P4 + k (4.6) 

we then have 

Z 2 e G m 2 M 2 v-^ fd s p 3 [<Pk, f ,, 2 . . 

da b = —, — ; > / — — / (5 (Yt — fe) — M ) #(e 4 ) 

(271") y/fa-piy-mW^J 63 7 W U ' 



,2 1 / P3 Pi „ P4 , „ P2 



x |«(p 3 )7^(pi)«(p4)r^(p 2 )r 7^2 M^-^T- + (4-7) 

(<?) VPs • k pi-k pA-k P2-k) 

in which p^ = p\ + p 2 — P3 — k. We may then transform to the special frame S° (defined by 
t = 0), in which the delta function in (|4.7]) is independent of the angle at which the photon 
is emitted. There 

(t - kf — M 2 = t 2 Q — 2t uj + A 2 - M 2 = 0; t = e x + e 2 - e 3 (4.8) 

The photon energy is then given solely by the final electron energy. The procedure used 
in flj is to integrate next over the photon energy and angle in S° and then transform back 
to the lab frame to integrate over the energy and angle of the final electron. Instead, we 
remain in the special frame and integrate first over €3, the delta function giving €3 in terms 
of uj. However, the range of photon energies is assumed to be sufficiently small compared 
to all other energies that we can set e 3 equal to its value for elastic scattering throughout 
the integrand. In addition, we take the angular range of the final electron to be sufficiently 
small that we can take some average value for these angles and neglect any variation of the 
integrand over this angular range. Similarly, we neglect k in the above expression for p 4 . 
We may then take \u(p 3 )'y fJ ,u(pi) u{pi)T^u(p2)\ 2 / (q 2 ) 2 outside of the integration, giving 

, Z 2 e e m 2 M 2 ^ 1 



2\2 



(2vr) 5 ^{ Pl .p2) 2 -m 2 M 2 (q 

x \ps[ r' fk (ji Pl z p 4 l z P2 y (49) 

2M J uj \p 3 ■ k pi ■ k Pa - k p 2 • k ) 



The term 2M in the denominator in (|4.9|) comes from the delta function in ( [4.7| ), which 
contributes the factor \d{5 ((t - k) 2 - M 2 )}/de 3 \ = 2\t - u\ from flOfi . A gain neglecting 
terms of order k, we note from ( [4.6[ ) that in S°, t — 00 = €4 = M. Comparing (|4.4|) , Q2.8| ) 
and (|2.10|) and noting that in arriving at ( |4.9|) we have neglected terms of order k in p 3 and 
p 4 , we have 

a [' d 3 k ( p 3 pi p 4 p 2 V 1n x 

= - — da / - Z + Z 4.10 

Air z J uj \ps ■ k pi • k p 4 • k p 2 • k J 



where uj = Vk 2 + A 2 . There then remains the integration over photon energy (restricted to 
|k| < Ae) and angles. The relevant integrals have been evaluated by 't Hooft and Veltman 
0, Sec. 7. We give here only their final result, rewritten using our metric; the essential 
steps in the derivation are given in their work. They define 
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d 3 k i 



in terms of which 



n , wrL n + m 2 L 33 - 2p x ■ p 3 L 13 

da b = -— da < +Z (-2p x • p 2 L 12 + 2p 3 ■ p 2 L 32 + 2pi ■ p 4 L u - 2p 3 ■ p 4 L u ) } (4.12) 



471-2 



+Z 2 (M 2 L 22 + M 2 L M - 2p 2 ■ p 4 L 24 ) 



As shown in ||, Sec. 7, for the case in which the momenta pi and pj are all on the mass shell, 
the integrals can, provided pi is not a multiple pj, be written in the form 

Ui = -, H' + gg'} (4.13) 



where 



and 



in which 



' Vi ■ Pj + yj (Pi ■ Pj) 



9 2 2 

A — mjm, 
1 j 



/2Ae\ 



s\f = 2 In | : — | In ( rr: ) (4.14) 



miirij J V A ) 



\miVt 2 / \rrijvt 2 
+L(1 _^ + Lfl m ' l - t 



^ Ti? / V Pi'fij 

( 3A-t\ ( m 2 l-t\ 

H'-y-H 1 ^) (4is) 



Pi -Pj + J (Pi - Pj) 2 - m 2 m 2 
a = 5 ; I = api — (4.16) 



A j = Pi j ■ * + \l (Pi,j ■ t) 2 - m 2 / 2 ; 



Hi 



\l (Pi ■ Pj] 



9 2 2 

z — mfm, 
1 j 



(4.17) 



The evaluation of ( 4.11 ) for pi = pj is straightforward. The result written in terms of 
relativistic invariants is 



47T 

2 

mf 



In 



2At 



Pi ■ t 



In 



\/ (pi ■ t) 2 — mft 2 yrriiy/t 2 



Pi 



(4.18) 



In 0, Sec. 7, the expression for S 1 ^ is evaluated in the frame S°, defined by t= 0. Since 
we want finally to express the cross section in terms of lab frame energies and momenta, we 
have, in ( |4.15| ), written Sff in terms of relativistic invariants. 
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The terms of leading order in In A are apparent in ( |4.14|) and (f4.18|) . Substituting these 
in ( [4. 12 ) gives the infrared divergent terms in do\>- They are cancelled exactly by the In A 
terms in the elastic cross section. 

We next express da^ in terms of lab frame energies. To that end, we assume that Ae 
is less than any of the other energies and therefore now neglect k in ( |4.6| ), taking p± to be 
given by its value for elastic (non-radiative) scattering: 



t = Pl +p 2 -p 3 =p 4 (4.19) 

(2) 

0, S^j remains finite; the only singularity is confined to the term 

7(1) 



(Note that for Ae 

ln(2Ae/A), evident in S^' .) The relativistic invariants in can then be written simply in 



terms of lab frame energies: 



Pi -P4 
P3-P4 

Pi -P3 
P2-Pi 



P2-P3 
P2-P1 



Me 3 

2 



(4.20) 



— q + m 
-\q 2 + M 2 



Further, we express Ae, the maximum momentum of the photon in the frame S°, in terms 
of the final electron detector acceptance in the lab frame, AE: 



Ae = r]AE 



(4.21) 



Details of the derivation are given in Appendix E. Substituting ( f4.19| ) and ( f4.20| j in ( [4.14[ ), 
( |4.15| ) and ([4.18|) then gives in terms of lab frame energies and momenta. We note that 



although Lij as defined in (|4.11| ) is clearly symmetric in % and j, the expression for in 



( |4.15| ) is not manifestly symmetric; the form of the expression for Ly is rather different from 
that for Lji. In writing the explicit expressions for the terms in L^, we choose i and j such 
that Lij simplifies readily for lab frame electron energies and momentum transfers which are 
very large compared to the electron rest mass. When irti ^ rrij, this is achieved by choosing 

(2) 

i and j such that rrii = m and nij = M (note below that we have written S^ 2 )• 

At this point we make the high energy approximation as defined in section II, in which 



case the above expressions for and S^' simplify considerably. For S^' this is straight- 



t(2) 



ij 



(2) 

forward. We give the results for Slj in Appendix F. Substituting these in L^, the high 



energy approximation for the inelastic cross section given in Q4.12 ) then becomes 

[in (=£)-! 
In 



da h 



a 



71 



In 



-<t 

m' 2 



- 1 



da < 



2 1/ 



+2Z 



£4 
|P 4 I 



\nr] In ( ^ 2n ff " > j — In rj \nx 
+L(1 - 2) - Ul - i) 

[ln 2 x - lnx + L(1 - 4 



(4.22) 
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V. RADIATIVE CORRECTIONS TO ELASTIC ELECTRON-PROTON 

SCATTERING 



The results given in ( |3.39|) , (|3.40|) , and ( |4.22|) may be added to give the radiative correc- 
tion, 8. The analytic expression is given below in Q5.1|) and (|5.2|) . Numerical evaluation of 
the radiative correction for various values of the pertinent parameters (initial beam energy, 
momentum transfer, final electron detector resolution, and target nucleus) are given in Ta- 
bles Q, O, and [TIT]. We note that the infrared (In A) terms, which appear in both the purely 
elastic ( |3.39| ) and inelastic fl4.22|) contributions to the radiative correction, cancel exactly 
when added to give the cross section for elastic electron-proton scattering with radiative 
corrections to first order in a: 



da = do"o(l + 8) 



(5.1) 



where 



8 = ^ 

71 



6 



+ 



+ 



2aZ 

71 

aZ 2 



71 



(1) 

el 



28 



In 



J^aeJ 



In ( 46163 



+ L 1- 



x 



L 1- 



1 

T]X 



- In rj 
2 ' 



R| ("I ln ' x - lux In + lux) - lnx - 1 j In 



1 



( M 2 



+ 1 



7l~ 



L(cos^ -6) - - 
(5.2) 



Here, 8^ is the contribution coming from the inclusion of form factors for the proton and 
integration over the entire range of four-momenta of the virtual photon in the proton vertex 
correction (see (|3.5| ), (|3.12|) , ( |3.13|) , (|3.25|) ); it is thus not included in the analysis given in 
HI and 0, denoted here as the soft photon approximation. Moreover, 8^ has no infrared 
divergent terms; these are all included in the soft photon approximation. 

In Tables |, 0, and |TTl| we compare the values of the radiative correction, 8, calculated 
in this paper (denoted by MTj) with those given by Mo and Tsai in for various kinemat- 
ics. The initial beam energies and momentum transfers have been chosen to correspond to 
experiments proposed or already performed at Jefferson Lab |]TTJ and SLAC ||12j| . The final 
electron detector acceptance, AE, has been taken throughout to be one percent of the final 
electron energy, €3. In the form factors (see( [2.5|) ), the parameter A has been chosen to be 
700 MeV/c throughout. The contribution of the terms in (O) are grouped according to the 
power of Z which appears there as a factor. The numerical value of each of these groups 
of terms is given in the rows denoted by Z ,^ 1 ,^ 2 . Values given in the column MTj in 
the row Z 2 do not include the contribution of the proton form factor (which are contained 
in 8^ ); they are given for comparison with the values in @. In the range of energies and 
momentum transfers considered here, the correction 8^' , due to the finite size of the nucleon 
(and integration over the entire range of four-momenta of the virtual photon in the proton 
vertex correction), is found in general to be much smaller than the other contributions with 
factor Z 2 , labeled explicitly in ( |5.2|) and in Tables [I], ||, and |HI| . The values given in these 
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tables include only the contribution of electron-positron pairs in the vacuum polarization; 
the contribution of muon and tau pairs is given by ( |A8| ) and QA10 ). 

The curves in Figs. |3| and |] illustrate the two aspects of the present work: (1) the 
contribution of nucleonic size effects to the radiative correction, and (2) the improvement 
of the mathematical treatment of the integrations given in the work of Mo and Tsai [|TJ, 
0. The nucleonic size effects are all contained in the term 8^\ (Eq. Q3.40 )); its contribution 
relative to the overall radiative correction factor, (1 + f^MTj), is given by the dashed curves 
marked VTX, where VTX = 100x8^/(1 + 5 MTj ). The dotted curve shows DO = 100x( 

^MTj — ^Tsai) /(l + #MTj)) which is that part of the difference between the radiative correction 
given by Mo and Tsai M and the one given in this paper due solely to the improvement 
of the mathematical treatment of the integrations. Here, is the radiative correction 

given in ( |5.2| ), excluding the term 8^'. It will be noticed that VTX is always positive, and 
that for most of the range of allowed momentum transfers, DO is negative. Thus their sum, 
which is the difference between the radiative correction ^MTj given in this paper in Q5.2D , 
and f^Tsai, given in ||, is rather small except for the region corresponding to large scattering 
angles. This sum is given by the solid curves marked D, where D = DO + VTX = 100 x( 

^MTj — ^Tsai)/ (1 + ^MTj ) • 



VI. CONCLUSION 

We have calculated the radiative correction to elastic electron-proton scattering to lowest 
order in a using a hadronic model which includes the finite size of the nucleon. The contri- 
bution from the emission of real soft photons by the electron and the proton is calculated 
exactly. The contributions of the box and crossed-box (two-photon exchange) diagrams are 
calculated in a soft photon approximation which is less drastic than that employed in 0. A 
number of observations may be made from the values given in Tables |, [H], and |TTT| . First, 
the contributions of the electron vertex correction, vacuum polarization, and real soft pho- 
ton emission by the electron (the terms in Q5.2| ) with factor a/n) dominate the radiative 
correction S. Since our expression for these terms differs from that given by Mo and Tsai 

solely in that they have omitted the term (a/n) L(cos 2 |(9) — ^ in ( |5.2| ), we find values 

for 8 which differ from theirs by at most 2% for the initial energies and momentum transfers 
considered here (note that — ^- < L(cos 2 |6>) — ^ < 0). Further, we note that, except for 

the proton and at the higher energies considered here, the contribution of 8^ is negligible. 
However, for the two highest energies, 5^ is between 2% and 3% of the factor (1 + 8) by 
which the uncorrected cross section must be multiplied, and hence should be considered 
in precision measurements for electron-proton scattering at energies above 8 GeV. As an 
empirical guide, we find that 8^ = 0.02(1 + 8) for initial energies and scattering angles 
satisfying e\ sin 6^8 for beam energies between 8 and 16 GeV. Finally, we note that a 
considerable simplification of the expression in Q5.2| ) occurs if, in addition to the last two 
terms multiplying a/n, we neglect the last two terms multiplying 2aZ/n as well as the last 
three terms multiplying aZ 2 /n, each of these sets of terms being always less than n 2 /6 in 
magnitude. From this study we see that at the energies and momentum transfers consid- 
ered here, the nucleonic finite size effects are rather small but are expected to become more 
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important at higher energies. The corrections due to the improvement of the high energy 
behavior of the radiative corrections as described in this paper are not negligible and need 
to be taken into account at the energies and momentum transfers we have considered. 
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APPENDIX A: ELASTIC SCATTERING AMPLITUDES TO ORDER a 2 

Using the notation given in section II, the matrix elements corresponding to the various 
one-loop diagrams shown in Fig. [l| are 

M 5 = Ze 2 u(p 3 )A,(p 3 ,p 1 )u(p 1 )^^u(p A )T^q 2 )u(p 2 ) (Al) 



M 6 = Z s e 2 u(p 3 )^u(p 1 )^-^-u(p 4 )A^p 4 ,p 2 )u(p 2 ) (A2) 

q \ 1£ 

where 

A M (p 3 ,Pi) = ~ie 2 [ 7^4 ] 2 , ■ Yji it 1 ^T7/,77 T~ T^l" ( A3 ) 

J (27r) 4 k 2 — X 1 + ie (p 3 — ft — m + ie) (pi — ft — m + ie) 



1 



M 2 = (Ze 



2\2 



d 4 k 



(fa - ft - M + ie 
1 1 



(j>4, - ft - M + ie 



-T^q 2 ) 



x 



(2tt) 4 k 2 - X 2 + ie (k + q) 2 - X 2 + ie 
u(P3)lv- n 1 - , ,, 7^(Pi) 



111 

pi — If — m + ie 



M 3 = (Ze 



2\2 



u(p^((k + q) 2 ) 
d 4 k 1 



1 



i> 2 + ft - M + ie 
1 



T»{k 2 )u{p 2 



(2tt) 4 k 2 - X 2 + ie (k + q) 2 - X 2 + ie 



u(pshu 



x 



u 



(p 4 )T»(k 2 



fx — ft — m + ie 
1 



^ 4 - ft - M + ie 



T»((k + q) 2 )u(p 2 ) 



(A4) 



(A5) 



(A6) 
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The matrix element, M4, for vacuum polarization is, after charge renormalization, related 
simply to the matrix element M\ given above: 



M 4 = U(q 2 )M 1 (A7) 
For a fermion loop in the photon propagator, IT(g 2 ) is given in || in terms of an integral 



which can be evaluated in closed form fFJ 



giving 



i% 2 ) = nV/m 2 ) - ^ { (1 - 1) VT+^iog 



VI + u + 1\ 5'. 

! 1 u - - > (A8) 



VTTu-lJ 3 



in which rrii is the mass of the fermion and u = 4m 2 /(— q 2 ). 
For —q 2 /m 2 ^> 1 this gives 



nV/™D = -|M-^r)-~!> (A9) 



7i [3 \mf J 9 

If we include the vacuum polarization amplitudes from particle-antiparticle loops of different 
masses, as has been done in several experimental analyses ||12||, then we have 



M A = M l J2 n/ {QVm 2 ) (A10) 

i 

In principle, once one includes particle-antiparticle pairs of mass greater than the electron 
mass in the vacuum polarization amplitudes, one should consider bosons as well as fermions. 
The matrix elements for vacuum polarization for a pair of structureless spin zero bosons in 
the closed loop, first given by Feynman ||13|| , may be found in a more accessible form in a 
paper of Tsai The result, corresponding to the equation above for fermions, is 

M boson = M 1 J2 n V/m 2 ) (All) 

i 

where 

nV«)^{i^io g (^±i)- M -i} (A 12 ) 



A more complete discussion of vacuum polarization should include a consideration of 
pion structure as well as the contribution of spin-one bosons, in particular the p meson. A 
detailed discussion of the hadronic contribution to vacuum polarization may be found in 
connection with calculations of the anomalous magnetic moment of the muon || and in 
connection with radiative corrections to high energy electron-positron collider experiments 

0- 
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APPENDIX B: PROTON VERTEX CORRECTION 



As noted in Sec. IIIA, the terms ggg,gsg, ... can be expressed in terms of the integrals 
given in (3.1). There, the integrals Iq, Jq, and K are scalars and hence are functions of the 
scalars p 2 ,, p|,and p 2 ■ Pi (and, of course, A 2 and A 2 ). Since we have on-shell particles in the 
initial and final states (p\ = p 2 = M 2 ), these integrals are functions of M 2 and q 2 . The 
integrals 1^ and are vectors and hence in principal can be written in the form 

(Bl) 



with a similar equation for J M , where a and b are functions of M 2 and q 2 . However, the 
calculation is simplified greatly if we express 1^ in terms of the four- vectors p = p 4 + p 2 
(which is symmetric in p 4 and p 2 ) and q = P4 — p 2 (which is antisymmetric in p 4 and p 2 ), 
i.e. 1^ = Ap^ + Bq^. Here A and B are functions of M 2 and q 2 and hence are symmetric in 
P4 and pi. Further, since the integrands for the vectors 1^ and are symmetric in p 4 and 
p 2 , it follows that B = 0. We thus have 



(B2) 



and a similar equation for J^. These same considerations of symmetry allow for the simpli- 
fication of the tensors and J M „, which are also symmetric functions of p 4 and p 2 . We can 
therefore write 



1 flV 



aip^pu + a 2 q^q u + a^g^ 



(B3) 



and a similar equation for J^ u . That the terms p^q v and q^p u are absent follows directly 
by multiplying successively by p^q v and q^p u , using p ■ q = and the fact that I fJjU p' 1 q lJ 
and I liv q il p v are antisymmetric in p 2 and p 4 . Multiplying J M (J^) by p M , and (J 1 ^) 
successively by p'V", g^g 1 ', and g^, the coefficients in the expressions for J M , J M , and J M! , 
may be expressed in terms of their moments, defined by 



0i 



p 2 



p 2 



(B4) 



0ii = VpV, 



fell = -^J^up^p" 



1 1 

022 = —^^(fcf, h 22 = -^J liV q t "q v 

Conversely, the integrals in ( p.l|) may be written in terms of the moments: 



Pn9i, 



Jfj, = Ppfrl 



(B5) 
(B6) 

(B7) 



IflU PfiPv 



0n ~y2^ ha ~ P 2 9u ~ ^922) 

022 - ^2 ~~ ^ n ~~ ^fe) 



:(^o - P 2 0n - <7 2 022) 



(B8) 
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J fjLU P^Pv 

+5V 



'ii 



2p- 



:{ko - p 2 h u - q 2 h 22 ) 



h 



22 



2q' 



■(k - p 2 hn - q 2 h 22 ) 



(B9) 



1 



(ho - p 2 h n - q 2 h 



22 i 



(BIO) 



where, for convenience of notation, we have defined 

9o = hi h = J , k = Kq 

The terms ggg,gsg, ... can be expressed in terms of these moments. Substituting ( |2.4| ) and 
( p.5|) in the expression for the proton vertex correction, Eq. QA4J) , and substituting this in 
turn in the matrix element M 6 , Eq. the integrals are all of the form given in Eq. ( |3.1|) , 

which are in turn expressed in terms of the moments by Eqs. ( p7|) -( pT0|) . We make use of 
the fact that in Mq the vertex correction is taken between free spinors, and finally express 
p 2 and P4 in terms of p and q. We find 



999 = -ie 2 F(q 2 ) 



2(2M 2 - q 2 )g - 4(2M 2 - q 2 ) 9l 
-2(8M 2 + q 2 ) 9ll - 2{q'/p 2 ) 922 + 8(M 2 /p 2 )h c 



7m 



+ [-8M 2 9l + 24M 2 9u + 8M 2 (q 2 / p 2 )g 22 - 8(M 2 /p 2 )h ] 



2M 



gsg = —ie 2 nF(q 2 



ggs + sgg = -2ie 2 nF(q 2 ) 



gss + ssg = —2ie 2 K 2 F(q 2 ) 



+ [2(2M 2 -g 2 )( 7 o-4(2M 2 -g 2 )^ 1 ] 



2M 



[-12M 2 g n - (q 4 /p 2 )g 22 + 2(1 + 2M 2 /p 2 )h - 3h x \ 7m 



+ 



(16M 2 - q 2 )g n + (g 2 /p 2 )(8M 2 - q 2 )g 22 
-4(1 + M 2 /p 2 )h + 3h 



-(8M 2 + q 2 )g n -(qyp 2 )g 22 
+ (-2 + 4M 2 /p 2 )h + h 1 
+ [3q 2 g n + AM 2 (q 2 / p 2 ) g 22 - (q 2 /p 2 )h - h x 



(q 2 /AM 2 ) 



2M 

7m 

2M 



(Bll) 



(B12) 



(B13) 



(B14) 



sgs = —ie k F(q 



-(8M 2 + q 2 )g n - (q 4 /p 2 )g 22 + (2 + q 2 /p 2 )h 
(8M 2 + q 2 )h u /4M 2 - (q 4 /p 2 )h 22 /4M 2 + (-1 + q 2 /p 2 )k /4M 2 
!2M 2 9ll + 4M 2 (q 2 /p 2 )g 22 - 2(1 + 2M 2 /p 2 )h 1 

+3h n +4(q 2 /p 2 )h 22 -k /p 2 ™ 



7m 



(g 2 /4M 2 ) 



sss = -ie K 6 F(q z ) < 



12M 2 9ll - 4M 2 (q 2 /p 2 )g 22 + 2(1 + 2M 2 /p 2 )h 

-Ah 1 + 3h ll + {q 2 /p 2 )h 22 -k /p 2 
2(2M 2 + q 2 )g n + AM 2 (q 2 /p 2 )g 22 
-(4M 2 /p 2 + q 2 /2M 2 )h + {q 2 lM 2 )h x 
-2(2M 2 + q 2 )h n /4:M 2 + 2{q 2 / p 2 ){2M 2 - q 2 )h 22 /AM 2 
+ (q 2 /p 2 )k /AM 2 



(B15) 

7m 



2M 



(B16) 
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The expressions do not depend on the particular form of the form factors; we have assumed 
only that F\ = F 2 = F. However, for form factors of the form given in Eq. ( |2.5|) , the 
moments may all be expressed more simply in terms of the functions C(A 2 ) : 

(C (A 2 );C M (A 2 );CV(A 2 )} = J d 4 k {1; k^ k^} /D(A 2 ) (B17) 

Using the identity 



k 2 - A 2 V^ 2 - A 2 / (m - 1)! 1 A 2 - A 2 
with T = Qfj^rj , we can write 



k 2 - A 2 k 2 - A 2 



(B18) 



{/} = N m (A 2 ) m T m ~ 1 { ° {A 2 _ gf^ } , (B19) 



where 



( —\) m 

i)W - (B20) 

In Eq. ( |B19| ) / and C denote any one of Iq,^, and C ,C^, C^ u respectively We see from 
Eq. (pi9|) that terms in C(A 2 ) which are independent of A 2 do not appear in the expression 
for / . In particular, this applies to C M ^(A 2 ), which may be evaluated using either dimensional 
regularization or a convergence factor. The infinities in C^ U (A 2 ) are indeed independent of 
A 2 , thus giving a finite result for as it should. In similar fashion, we have 

{J} = AC(A 2 rT— { A2C(A A 2 2 } :^ (A2) } (B21) 

in which J and C denote any one of Jo, J M , J^u and Co, C M , C^ u respectively. We see from Eq. 
( B21|) that any terms in C(A 2 ) which are independent of A 2 do not appear in the expression 
for J provided m > 1. Finally, for K we get 

Ko = N' m ^TT^ { A 4 go(A^ 2 )-An(A 2 ) j ^ 



We note that, apart from trivial factors, the integrals in (|B17| ) are the three-point func- 
tions defined in ||, Eq.(5.1), and ||, Eq.(E.l); C has been evaluated in terms of Spence 
functions in ||. The details of the algebra in || and [|J being rather lengthy, we choose 
instead to evaluate the integrals in (P17| ) using Feynman parameters, writing 

1 ■>!"[:. ... (B23) 



D(A 2 ) J J [k 2 - 2xk ■ Py - A 2 (l - x) + ie] 

where p y = P2V+P^ — y)- Substituting Eq. ( B23|) in Eq. (|B17|) and shifting the integration 
variable k (k — xp y — > k) we then have 



21 



|x, a; Py/ti 2< PyiiPyu 4" n^ui/A {x ) | 
^d?/- ^9 , A 9/, 3 ^ + 



-27T 



^0 



x 2 p 2 + A 2 (l - x) 



(B24) 



where, throughout this section, we denote by x^v any terms which are independent of A . 
Rewriting p y = \p + \q(l - 2y) in Eqs. (|B23| ) and ( ggg ) we may, neglecting terms which 
are independent of A, express Co, C M and in terms of the functions 



We get 



1 /■! 



cix 



./0 



p 2 o; 2 + A 2 (l — x) 



C o = -m 2 0i(A 2 ); C M = -m 2 -p^ 2 (\ 2 ) 



(B25) 



(B26) 



C 



-27T 



1 lp 2 



A 2 [<MA 2 ) - 2 (A 2 )] + i^A 2 [0!(A 2 ) - ^0 2 (A 2 )] 



(B27) 



As shown in Appendix D, the functions <pk obey a three-term inhomogeneous recursion, 
which is used to calculate fa for k > 1: 



(& + l)p 2 fc+2 (A 2 ) - 2(2A; + l)A 2 fc+1 (A 2 ) + 4£A 2 fe (A 2 ) 

ki(A 2 : 



pi 



(P + Pl\ 


+2A 2 


\p-pl) 





Here 



^ o) (A 2 )=0 fe (A 2 )i_ 92=o = r 

^0 



k + 1 



/(! M 2 a; 2 + (1 - a;)A 2 
The functions ^^(A 2 ) may in turn be calculated from the recursion 

M 2 ^(A 2 )-A 2 4°] 1 (A 2 )+A 2 0f(A 2 ) 1 

To implement the recursions ( |B28| ) and ( |B30| ) we need 

4 0) (a 2 : 

and 

#(A 2 ^- 1 



1 , fA + A x 
in 



AAi V A - A i 



2M 2 



' M 2 A /A + Ai 

ln it + i - ln i r~ 

A 2 Ai VA-Ai 



(B28) 



} (A 2 



(B29) 



(B30) 



(B31) 



(B32) 
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which follow from Eq. (P29|) , and 0i(A 2 ), which can be expressed in terms of dilogarithms 
(Spence functions) (see Appendix D): 

0i(A 2 ) = — [l (l — ^ -21n(x)ln(l + i)l (B33) 

ppi I V xv J V yJ y J 

where 

L[z)= _ n<LlA dt (B34) 



X= £±£L =<£+£>! (B35) 

P — P\ 4:M Z 

A + Ax (A + Ai) 2 

^a^aT^m^ (B36) 

We will also want to take the limit A —>■ 0. Neglecting all terms which vanish in this limit, 
we find, 

^ . 1 f or ( X \ ^ , f P 2 



and for k > 1, 



A -> OM 2 V A 



■(") = (T^r 1 " (B39) 



4 ™<°» = (F^ (B«) 



We then find 



#o = -iVi0i(A 2 ) + N m (A 2 ) m T m ~ 1 j-l^A 2 )} (B41) 
9l = ^N m (A 2 ) m T m ~ 1 |^ [0 2 (A 2 ) - MO)] } (B42) 
0u = ^Am(A 2 ) m T m_1 j-i [0 3 (A 2 ) - 3 (O)] } (B43) 



+^- 2 N m (A 2 ) m T m - 1 |0i(A 2 ) - l -H^ 2 ) 
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922 = -^N m (A 2 ) m T m ~ 1 ' 1 



A 2 



[0 3 (A 2 ) - 03(0)] 



-^AUA 2 )"^- 1 j<MA 2 ) - ^0 2 (A 2 ) 



(B44) 



/ io = iV m (A 2 ) m T m - 1 {0 1 (A 2 )} 



hi — -A^ m (A 2 ) m T m_1 (0 2 (A 2 )} 



ft-n = -AT m (A 2 ) m T m_1 |0 3 (A 2 )} 



+— N m (\ 2 ) m T m ~ l { A 2 



^(A 2 ) - -0 2 (A 2 ) 



(B45) 
(B46) 

(B47) 



h 2 2 = -f^N m (A 2 rT m - 1 (0 3 (A 2 )} 



-— ^(A 2 )^- 1 <| A 2 



> X (A 2 ) - jW) 



(B48) 



where 



A; o = iV m (A 2 rT m - 1 {A 2 1 (A 2 )} 



(B49) 



(B50) 



The terms ggg,gsg, ...sss can now be expressed more simply in terms of the functions (pk- 
We get 

-2(2M 2 -g 2 )AT 1 1 (A 2 ) 

ggg = -ie 2 F{q 2 ) { +2(2M 2 - q 2 ) [S^ 1 {^M 2 )} ~ S m ' 1 {± [0 2 (A 2 ) - 2 (O)]}] } ^ (B51) 

{0 2 (A 2 )} - 4M 2 S"- 1 [0 3 (A 2 ) - 3 (O)]} 



*e 2 F(g 2 ) \ —4M 2 S m ~ 1 \ - 2 [0 2 (A 2 ) - 2 (O)] \ + AM 2 S m ~ l \ - 2 [0 3 (A 2 ) - 3 (O)] 



A 2 



2M 



gsg = - ie 2 KF(q 2 ) j-g 2 ^ 1 j ^ [0 2 (A 2 ) - 2 (O)] } } 7 , 



(B52) 



-ie 2 nF(q 2 ) |2(2M 2 - g 2 ) 



-jV^A^ + ^j^MA 2 )} 



2M 



™ S + S ™ = -2^F(, 2 ) J & ~ 3M 2 )S"^ [0 3 (A 2 ) - 3 (O)]} 
^ + ^ 6 + |5— 1 {0 1 (A 2 )-i0 2 (A 2 )} ^ 

2M 



(B53) 



+ 



- Ke\F(a 2 ) I 2M2Sm ' 1 [03(A2) ~ ^ 
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+ = -2,eVF(g 2 ) — | _^ m _ 1 j 7, (B54) 

+ 2re 2 n 2 F(a 2 ) I - M2)Sm ~ 1 & [03(A2) ^ 1 ^ 

f — 2M 2 S ,m_1 [0 3 (A 2 ) - 3 (O)]} - {0 3 (A 2 )} ) 

sgs = -ie 2 K 2 F(q 2 ) { +S m ^ {0i(A 2 )} + \S m - x {0 2 (A 2 )} ) ^ (B55) 

I -^- 1 {A 2 [0i(A 2 )-|0 2 (A 2 )]} J 

2M 2 ^- 1 {^[0 3 (A 2 )-0 3 (O)]} li^9" 
6 [q } { -2S m - 1 {0!(A 2 )} + iS™- 1 {0 3 (A 2 )} J 2M 

. 23ri/a v g 2 f -2M 2 S m - 1 {-^[0 3 (A 2 )-0 3 (O)]} + iS m - 1 {0 3 (A 2 )} 1 ,_, ceN 

r ^^{^[^(A^-^wn-a+^^^^A 2 )} i „ 



-te 2 K 3 F(q 2 ){ -i5 m - 1 {^(A 2 )} 



■1(1 + A )^ 1 {^(A 2 ')} - {A 2 [0i (A 2 ) - 2 (A 2 )]} 



2M 



where 



= N m (A ) m T m ~ (B57) 



It should be noted that the terms with 0i(A 2 ), which appear only in ggg and gsg, 
constitute the well-known infrared divergence. They are, apart from the hard-photon proton 
interaction, ( [2.4|) , independent of the proton form factor (in this case independent of A and 
M). This is to be expected, since this term is cancelled by a similar infrared divergent term 
coming from the cross section for the emission of a real soft photon, which is given by the 
elastic cross section multiplied by a factor independent of the proton form factor. 



APPENDIX C: ELECTRON VERTEX CORRECTION 

As noted in Sec. IIIC, the electron vertex correction, M5, may be obtained directly from 
the proton vertex correction, M 6 , if we retain only the term ggg, set F = 1, replace p 2 ,P4 
and M by pi,p 3 and m and take the limit A — > 00. After making these replacements in 
ggg as given in ( P51|) , we need 0fc(A 2 ) to order m 2 /A 2 and </>fc(0) (see (|B39|)). From ( |B33| ) . 
writing ( §, p. 389 (B.2)) 

L(l - z) = -Liz) + -7T 2 - In* ln(l - z) 

6 

and neglecting terms of relative order m 2 /A 2 , we have 



and from (|B25|) 
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fc (A 2 )-0 fe+1 (A 2 



1 f 1 



A — >■ oo A 2 /c + 1 



Choosing m — 1 in (|B51|) we then have 



000 



2m 



where 



G (e) (g2) = « I _ 2(2m 2 _ g 2 )0i(A 2 } + [3pl 4m 2 



47T 



PmPl 



1 , f A " 

In x m + - + In — 

2 V mr 



(C2) 



(C3) 



(C4) 



a f 4m 2 



In x r . 



47T tp m pi 

Adding the contribution of the electron self energy diagrams gives 

m = [(Gi e V) - Gi e) (0)) 7|t + 4V)^f 



(C5) 



(C6) 



To facilitate comparison with [T|], we write this in terms of the functions K(pi,Pj). Similarly 
to (13.71), we now have 



(C7) 



which then gives ( 3.18 ). 



APPENDIX D: THE FUNCTIONS fc (A 2 ) 

In this Appendix we derive the three-term recurrence relation for the function 0^(A 2 ) 
given m ( |B28D as well as the expression for the function 0i(A 2 ), defined in ( |B25|) and given 
in terms of Spence functions in ( |B33| ). Integrating first over y in Eq. ( |B25| ) we have 



0,(A 2 ) = - / V lo § 
Pi Jo R 



R + xpi 
R — xpi 



(Dl) 



where R 2 = p 2 x 2 + 4(1 — :r)A 2 and p 2 = — q 2 > 0. Noting that 

4- \x k R\ = x k ~ l {kR 2 + p 2 x 2 - 2xA 2 } R~\ (D2) 
ax 1 J k J 

we get 

(k + l)p 2 <P k+2 - 2{2k + l)A 2 fe+1 + 4M 2 fc = - / log S^PlX d . ( D3 ) 

Pi Jo I K ~ x Pi ) 



Integration by parts then gives 
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(k + l)p 2 <p k+2 - 2(2k + l)A 2 fc+1 + 4A;A 2 fc = log 



2p, /p + Pi 



Pi 



2A 2 



x fc (2 - x) 



o M 2 x 2 + (l-x)A 2 
(D4) 



from which the recursion ([B28 ) follows at once, using Eq. ( |B29|) 

From Eq. ( p4|) it is clear that cf>2 and 03 follow once we have evaluated 0i, which follows. 
Setting = 1 in ( |B25| ) and integrating first over x, we get 



x dx 



o p 2 x 2 + A 2 (l-x) 2p 2 y 



2 

A 2 



h i I ft 1 +^ln 





4)] 







(D5) 



where A 2 = A 2 — 4p 2 . We next make the change of variable y — (1 + ^)/2, which gives 
A 2 = p\u 2 + A 2 — p 2 , and then make the further change of variable A = piuj + s, from which 



to 



A 2 - p 2 - s 2 
2pis 



A 2 - p 2 + S 2 
2s 



(D6) 



Then integrating ( |B25|) over y gives 



0i (A 2 



In 



PiV s _ P 2 -(A-s) 2 
2 r+ ds 
pT7 s _ (A + s) 2 - p 2 



(A + s) 2 - p 2 



In 



4sA 

P 2 -(A-^) 
4sA 



(D7) 



where 



s± = Ai ± pi 



;ds) 



Factoring the expressions which appear as factors to the logarithms as well as in their 
arguments, we can write 



where 



h 



ds 



s — er_ 



>i(A 2 ) = — x> 



■ln(s + cr + ); I 2 



ds 



S + O- 



ln(s — cr_ 



(D9) 



ds 



OA 



ln(s + cr_); I4 



(is 



S + <7_ 



In (a+ 



ds , / s + cr_ 

m I — | ; i 8 



s — <r_ 



2s 



(is / -S — (T_ 

m 



S + cr_ 



2s 



(D10) 
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s + a i 



2s 



G?S . ( S + CT-i 

In 



2s 



ig 



1 



1 



1 



1 



ln(2A) 



s — er_ o~ + — s s + cr_ s + cr H 
where 

cx± = A±p (Dll) 
Some of the integrals J, can be integrated out directly. For these terms we get: 

4(74. cr_ 



c4)(l-ai) 



lna+ln(l-a;)+lna_ln(l-ai) (D12) 



and 



J 9 = -2 In — ln(2A) 



(D13) 



where 



a± 



PTPi 



A + Ai 

The remaining terms in Eq. (P9|) we may rewrite as 



In = In «+ ln(l - c^) - In a_ ln(l - a 2 _ ) + L(l - a^) - L(l - a 2 _ ) (D14) 



n=5 



where L is the dilogarithm (Spence) function. We then get the result given in ( |B33| ) in 
Appendix B. 



APPENDIX E: FINAL ELECTRON DETECTOR ACCEPTANCE 

In this Appendix, we express Ae, the maximum momentum of the photon in the frame 
S°, in terms of the final electron detector acceptance in the lab frame, AE. In S° (p4+ k 
= 0), if |k| = Ae M, we have, from (pi + P2 — P3) 2 = (p4 + k) 2 , neglecting terms of order 
(Ae/M) 2 and (m/M) 2 , 



P2 ■ (pi - P3) ~ Pi ■ P3 = MAe 
Writing this in terms of lab frame energies, we have, for high energies, 

M(e x - e 3 ) - 6x63(1 - cos#) = MAe 
For elastic scattering in the lab frame, we have 



(El) 



(E2) 
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M(ei -ef) -eiegCl -cosfl) = (E3) 

Subtracting gives 

A£ (l + -^-(1 - cosfl)) = Ae (E4) 



M 
where 

AE = ef - e 3 (E5) 
Thus, in terms of lab frame quantities we have 

Ae = 7]AE (E6) 



APPENDIX F: HIGH ENERGY APPROXIMATION FOR S, 



(2) 



(2) 

In this appendix we give the high energy approximation of the terms defined in 
( |4.15 ), in which we note in particular that for % = 1 or 3 we have l-t = (api — pj) -t ~ api ■ t. 



Using transformations of the dilog (Spence) functions [[|, p. 389 (B.3), 

L{z) = -l(-\ - l -\n\-z) (Fl) 



L{z) = -L (-iL-J - iln 2 (1 - z) (F2) 

(2) 

the terms in Sh simplify considerably. We then obtain 



S[f = - In 2 l 2 ^) - In 2 x + J In 2 f ^ - ivr 2 (F3) 



TO / 2 \r] 

xr] / V xJ 



S® = - In 2 f^j - In 2 (*) + ~ In 2 (x V ) - (F4) 



+Hi-1-)-l(i-V 



XT] J \ X 



Sgl-.tf » (F5) 

TO / O 
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s 



(2) 
13 



In 2 



2ei 
m 



2 ^2e 3 



: / i cos 2 -6 



In 2 

\ m 



1 2 1,2 
7T 2 + - In 2 

3 2 



cos 



- 1 * 



(F7) 



S. 



(2) 
2d 



In 2 (aO 



-L I 



(F8) 
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FIG. 1. Feynman diagrams for elastic amplitudes. 



FIG. 2. Feynman diagrams for inelastic amplitudes. 



FIG. 3. The curves show the contribution of nucleonic size effects (VTX - dashed curve), math- 
ematical refinement (DO - dotted curve), and the resulting difference (D - solid curve) between the 
radiative correction given by Mo and Tsai [Q] and that given in this paper, 5mtj j as a function of 
four-momentum transfer for an initial electtron energy of 6 GeV. (VTX, DO, and D are defined in 
Sec. V.) 
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FIG. 4. As Fig. H|, but with an initial electron energy of 16 GeV. 

TABLE I. Contributions to the radiative correction 5 for electron-proton scattering as given 
in this paper (MTj) and in Mo and Tsai Q (MoTsai) for three initial electron energies and 
four-momentum transfers. Values in the rows marked Z ^ 1 , and Z 2 refer to contributions from 
terms with these factors in (|5.2|). 



ei = 4.4 GeV 




ei = 12 GeV 




£l = 21.5 GeV 




Q2 = 6 (GeV/c) 2 




Q 2 = 16 (GeV/c) 2 


Q 2 = 31.3 (GeV/c) 2 


MTj 


MoTsai 


MTj 


MoTsai 


MTj 


MoTsai 


Z° -0.2187 


-0.2171 


-0.2330 


-0.2322 


-0.2323 


-0.2317 


Z 1 -0.0569 


-0.0506 


-0.0517 


-0.0479 


-0.0625 


-0.0571 


Z 2 -0.0242 


-0.0232 


-0.0359 


-0.0347 


-0.0452 


-0.0440 


$ +0.0068 






+0.0116 




+0.0185 




5 -0.2930 


-0.2908 


-0.3090 


-0.3149 


-0.3214 


-0.3328 


TABLE II. Contributions to the radiative correction 5 as given in this paper (MTj) and in Mo 


and Tsai M (MoTsai) for several nuclei, with e± = 4.4 GeV, Q 2 


= 6 (GeV/c) 2 ; other symbols as 


in Table |[ 














2 H 




4 He 


12 C 


40 Ca 


MTj 


MoTsai 


MTj 


MoTsai 


MTj 


MoTsai MTj 


MoTsai 


Z° -0.2476 


-0.2467 


-0.2535 


-0.2532 


-0.2615 


-0.2609 -0.2632 


-0.2625 


Z 1 -0.0187 


-0.0183 


-0.0066 


-0.0077 


-0.0151 


-0.0168 -0.0147 


-0.0173 


Z 2 -0.0094 


-0.0088 


-0.0077 


-0.0071 


-0.0156 


-0.0145 -0.0188 


-0.0178 


$ +0.0010 




+0.0002 




+0.0001 


+0.0001 




5 -0.2747 


-0.2739 


-0.2677 


-0.2680 


-0.2920 


-0.2922 -0.2966 


-0.2975 



TABLE III. Contributions to the radiative correction 5 as given in this paper (MTj) and in Mo 
and Tsai ^ (MoTsai) for several nuclei, with e\ = 21.5 GeV, Q 2 = 31.3 (GeV/c) 2 ; other symbols 

as in Table | . 

2 H 4 ~He ™C 4^ 





MTj 


MoTsai 


MTj 


MoTsai 


MTj 


MoTsai 


MTj 


MoTsai 


z° 


-0.2707 


-0.2704 


-0.2817 


-0.2814 


-0.2876 


-0.2875 


-0.2896 


-0.2894 


z 1 


-0.0182 


-0.0189 


-0.0149 


-0.0166 


-0.0134 


-0.0164 


-0.0123 


-0.0173 


z 2 


-0.0231 


-0.0221 


-0.0379 


-0.0352 


-0.0583 


-0.0535 


-0.0761 


-0.0708 


® 


+0.0045 




+0.0026 




+0.0006 




+0.0001 




s 


-0.3076 


-0.3114 


-0.3319 


-0.3333 


-0.3587 


-0.3573 


-0.3786 


-0.3775 
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elastic electron electron vertex electron self-energy vacuum 

scattering correction diagrams polarization 




proton vertex proton self-energy box and crossed- 

correction diagrams box diagrams 
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